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- - - Abstract 

In this paper the longitudinal profile of muon production along the shower axis is studied. The 
^ I characteristics of this distribution is investigated for different primary masses, zenith angles, pri- 
\ mary energies, and different high energy hadronic interaction models. It is found that the shape of 
■ this distribution displays universal features similarly to what is known for the electromagnetic pro- 
file. The relation between the muon production distribution and the longitudinal electromagnetic 
OO . evolution is also discussed. 

■ . Keywords: Extensive Air Shower, Longitudinal Profile, Muon Production Depth, 
p I Electromagnetic component 

1. Introduction 

^ ! The origin and mass composition of the most energetic particles in the Universe, the Ultra 

Tj- ! High Energy Cosmic Rays (UHECRs), remains a mystery. These particles reach the Earth with 
^ \ a very scarce flux. Fortunately, their interaction with the atmosphere molecules produces huge 
cascades of particles, known in the literature as Extensive Air Showers (EAS). The detection of 
these showers is usually done by measuring the charged particles that arrive at the ground, or, in 
moonless nights, the development of the shower can be followed through the fluorescence light 
produced by the EAS. 

While the arrival direction of the UHECRs can be easily obtained using any of the techniques 
mentioned, the mass composition of the primary particle is much more difficult. The shower 
observables connected to the type of primary particle are also sensitive to the physical interactions 
5^ . that occur during shower development. The hadronic interactions at high energies are described 
through phenomenological models that are fitted to the available accelerator data and extrapolated 
several orders of magnitude to the UHECRs energies. Moreover, the accelerators have difficulty 
to reach the most important region for the characterization of the EAS development, the forward 
region, increasing the uncertainties on the extrapolation. 

Contrary to the fluorescence light, which is dominantly produced by low energy electrons in 
secondary electromagnetic cascades, ground signals are sensitive to muons produced at different 
depths, thus imaging the hadronic cascade. 
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Many of the muons decay before reaching the ground. Current EAS Monte Carlo simulations 
(for instance CORSIKA ^ or AIRES [2]) properly account for the propagation effects starting 
from the moment of production. The total/true Muon Production Depth distribution (MPD) is de- 
fined as the total number of muons produced in each slant depth unit, regardless of the probability 
to reach ground and be detected. On the other hand, the apparent-MPD distribution is affected 
by the geometrical and propagation effects and it has been extensively studied aaaaai. 
In 0] the apparent is related to the true MPD distribution through the energy, transverse momen- 
tum at production and propagation effects. This knowlegde would allow the reconstruction or at 
least constrain the total/true MPD, which is closer to the development of the hadronic cascade, i.e. 
without being masked by propagation effects. 

The absolute density of muons at ground has been used previously to study high energy 
hadronic interaction models and primary composition of cosmic rays. The number of muons 
at a given distance from the core is analised in [10], showing that at very high energy the data 
can not be explained by the available models. The simultaneous analysis of the electromagnetic 
calorimetric energy and ground signals has been used in [11] to correct for the invisible energy 
carried by neutrinos (and muons), showing that the energy determination can be improved on an 
event-by-event basis. 

In this paper we identify the main characteristics of the total/true MPD distribution charac- 
teristics at UHECRs energies. The shape of the depth profile, depth of maximum, X^ax^ and its 
corresponding value, N^ax will be studied. The comparison with the features observed in the en- 
ergy deposit profile, referred to in this paper as electromagnetic profile, will be made whenever 
relevant. 

Since in this study we are only interested on the longitudinal profile, and not the transverse 
distributions, we used the hybrid shower simulation CONEX Il3n (version v2r3). This program 
combines Monte Carlo simulation with one-dimensional cascade equations making it very fast and 
allowing the production of large samples of showers. CONEX gives as output the energy deposit 
profile as a function of depth (X), the number of muons produced as a function of X, and the more 
extensively used number of muons along the shower axis. Note that this last one, the number of 
muons as a function of the shower depth, is in first approximation the cumulative of the number 
of muons produced, minus muons that decay. Instead we will concentrate our study in the muon 
production profile. 

In this work we study the dependence of the muon shower production profile on: the primary 
mass composition (proton and iron) at several energies; the zenith angle at £■ = 10^^ eV; different 
hadronic interaction models, in particular QGSJet-II.03 [14, 15] and EPOS 1.99 fl^ at fixed energy 
E = 10'^; the primary energy (from log(£'/eV) = 17.5 - 20.0 in steps of 0.5). For each different 
set of parameters samples of 50 000 showers were generated. The high energy hadronic interaction 
model used as default was QGSJet-II.03. The ground was set at 4000 gcm^ to avoid an abrupt 
termination allowing us to see the full tail of the shower profile as if it had a very inclined zenith 
angle, 6. 

The paper is organized as follows: in Sec. |2]the main features of the muon production shower 
profile are studied and related to the electromagnetic profile; the shape of the profile is character- 
ized in Sec. Eland here it is shown that, similarly to what happens in the electromagnetic case, 
this profile exhibits an universal behavior when expressed in the coordinates, X' = X - X^ax and 
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A^' = NINmax', in section|4]the impact of the muon energy threshold on the longitudinal profile is in- 
vestigated; the extraction of more information on basic variables like the total number of muons or 
the point of first interaction are discussed in Sec. [51 the paper ends with conclusions and prospects. 

2. Longitudinal shower profiles 




Figure 1 : Muon production shower profiles as a function of depth (X), for proton (red) and iron (blue) primaries at 
E - 10'^ eV. In this picture is shown 100 showers for each primary, both generated with QGSJet-II.03 and with 
= 40°. 



The muon production longitudinal profile for proton induced showers (in red) and iron pri- 
maries (blue) at £■ = 10^*^ eV is shown in figure [T] These profiles were obtained with CONEX 
where the minimum energy threshold for muons is 1 GeV (the effect of this cut is discussed in 
Sec. SI). 

The shower profiles reflect the properties of the first stages of the hadronic interaction, in 
particular the first one, which can not be observed directly. They are thus very different for proton 
and iron initiated showers, as shown in Fig. [T] Compared to proton showers, iron showers have 
a higher number of muons, which is readily seen just by looking at the maximum of the profileQ, 
N^nax^^ iron showers more charged pions are produced and as a consequence it has more muons. 
Low energy pions will give rise to muons as they decay, while the pions above the critical energjH 



in this paper the index ji and e.m. will be used to address to variables related with the muon production profile 
and the electromagnetic profile, respectively. 

■^the critical energy is defined as the energy where the decay lenght equals the interaction length. 
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Figure 2: X^a c distributions of muon production longitudinal profiles for: (a) different primaries - the proton distribu- 
tion is in red/full while iron is in blue/dashed (showers generated with QGSJet-II); (b) high energy hadronic interaction 
models - QGSJet-II is the red/full line and EPOS 1.99 is shown as blue/dashed (for proton induced showers). 



are more likely to interact producing more pions. Protons exhibit much more fluctuations: this 
is due to the smaller pion multiplicity in the first interactions, thus less muons are produced and 
these are more spread along the shower axis. 

Neutral pions of all energies decay immediately and around 90% of the primary energy will 
feed electromagnetic cascades; the energy deposit in the atmosphere is dominated by low energy 
electrons after a few radiation lengths. All of the sub-cascades add up to give a well defined maxi- 
mum height of the profile. Not only the integral of the profile is a good calorimetric measurement 
of the energy of the primary particle (which must be corrected for the muons and accompanying 
neutrinos which reach the ground), but also the maximum of the profile is directly proportional to 
the energy, within 5% [17] for both primaries. Because it corresponds to a much higher number 
of particles than those in the muon production profile this number has much less fluctuations. In 
the electromagnetic case, it is the depth of the maximum, X^"'^, that gives more information about 
the other shower properties. 

The depth of the maximum of a shower, X,„ax, is determined by the depth of first interaction, 
Xi, but also the subsequent development, AX. The variation of the first is common to both the 
electromagnetic and muon profiles, while the later is diff"erent. L^on primaries have a larger cross- 
section and higher multiplicity and both eff"ects contribute to make the average maximum depth, 
(Xmax), and RMS (X,„ax) smaller than the corresponding proton values. In Fig. |2l we show the 
X'^ax distributions for the muon production profiles, comparing difl'erent primaries and difl"erent 
hadronic interaction models. The features are similar to the ones of the electromagnetic and 
the correlations between muonic and electromagnetic maxima can be seen in Fig. [31 together with 
the difl'erences. 

The electromagnetic maximum is reached around 200 g cm~^ later, after the energy is degraded 
to a large number of electrons (which will then stop multiplying), while the muonic profile has a 
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Figure3: Relation between muon production X^ax and electromagnetics^™^, for £ - lO'^eV showers. The difference 
between the two X,„av are shown in the inset plots. In (a) the dependence on primary mass is shown (proton is the 
red (full) line and iron the blue (dashed) line, both generated with QGSJet-II), while in (b) the dependence on high 
energy hadronic interaction models for proton induced showers (QGSJet-II is the red (full) line and EPOS 1 .99 the 
blue (dashed) line). 



maximum when the hadronic shower is still important. There is a correlation between the elec- 
tromagnetic and the muonic X,nax, shown in fig. |3l It depends slightly on the primary mass and is 
almost independent of the hadronic interaction models. 

The difference between electromagnetic and muonic X,„ax, shown as an inset plot of fig. |3l 
also changes between primaries - by around 30 g cm~^ in the average separation - which means 
that although related, the electromagnetic and muonic profiles give some independent information. 
Notice that this difference, X,^,™^ - X^qjc, is only sensitive to AX, thus it is giving information about 
the history of pion production. 

3. Universality of the longitudinal profile shape 

In Fig. |4l the profiles are expressed in X' = X - X,„ax and A'^' = N/Nmax- The obtained shape 
is rather universal, similarly to what happens to the energy deposit profile. It can be useful to use 
the average shape in order to determine the two main parameters from a fit with a reduced set of 
data. On the other hand, there can be extra information on the shape, and we can now construct 
average profiles to look in detail for differences between primaries and hadronic interaction mod- 
els. From that comparison, in fig. [5] (left), it is clear that iron showers develop faster, with almost 
no difference between models. 

The dependence on zenith angle is also studied. Muon production depends on the competition 
between the pion interaction and decay lengths. The larger the energy the most probable it is for 
a pion to interact instead of decaying into a muon, but that dependence can not be expressed in 
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Figure 4: Muon production shower profiles from proton (red) and iron (blue) primaries, in (X', N') coordinates. The 
same showers used to build Fig. [T]are used here. 



X alone, as the decay length is independent of the atmospheric density. Nevertheless, the average 
muon profile, in figure [5] (right), is still rather universal. There are, as expected, some minor 
violations to this universality that are more important in the early phase of the shower. 

This quasi-universal shape is compared to the electromagnetic shape in fig. [6l The diff"erences 
are clear: the muonic profile has a steeper growth and is more asymmetric, with respect to the 
shower maximum. Both profiles are fitted with a Gaisser-Hillas function, written in terms of 
a Gaussian width L and a parameter, R, which is related with the asymmetry of the shower in 
respect to the shower maximum iITtI] . 



We conclude that the muon production profile can be well described with the same function as 
used for the electromagnetic profile, as seen in figured It should be noted that the region around 
the maximum is better described by a Gaisser-Hillas function for the electromagnetic profile. How- 
ever, the full profile description is better achieved through the Gaisser-Hillas parametrization for 
the muon production profile. This is because the end tail of the energy deposit profile has impor- 
tant indirect contributions from muon decays. For showers initiated by protons of 10^^ eV, the 
width of the average muon production profile is larger by 40 g cm"^ and the asymmetry almost the 
doubled with respect to the electromagnetic profile. The average values and dispersion of these 
parameters evolve slowly with log(£'), as shown in fig. |7l for both primaries (proton and iron). 




(1) 
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Figure 5: Average muon production longitudinal profile in (X',N') coordinates: (left) dependence on primary mass 
- proton (red) and iron (blue), showers generated with QGSJET-II at 6* = 40°; (middle) dependence on the hadronic 
interaction model - QGSJet-II (red) and EPOS 1.99 (blue), for proton induced showers atO = 40°; (right) dependence 
on the zenith angle, - Black for 0° < 6i < 10°, red for 30° < < 40°, and blue for 45° < 6* < 55°, for proton showers 
using QGSJet-II. 
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Figure 6: Average shower profiles for proton primaries at E = 10'^ eV, with QGSJet-II, in (X',N') coordinates. 
Comparison between electromagnetic (in red) and muonic (in blue) shape features. The lines correspond to fits using 
a Gaisser-Hillas function (2 parameters). The fit results are given in the plot for the electromagnetic (e.m.) and muonic 
(fi) profiles. 
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The differences found for proton and iron profile shapes can now be quantified in these param- 
eters, with the corresponding distributions shown in fig. [81 The asymmetry is almost the same for 
both primaries, the diff'erence in means being well below the 5% dispersion in each sample. Since 

only aff"ects the tails it can easily be fixed in the following analysis. The Gaussian width, on 
the other hand, has diff"erent means for each primary, well above the single primary dispersion, 
and consistently for the two hadronic interaction models studied. So, it is a new variable for mass 
composition studies, and fairly model independent. 




17.5 18 18.5 19 19.5 20 17.5 18 18.5 19 19.5 20 

log(E/eV) log(E/eV) 

(a) (b) 

Figure 7: Shape parameters dependence on the shower energy. In (a) is shown the L parameter while in (b) is the results 
for R. The shape parameters for the electromagnetic profile are shown in red (full) line while the muon production is 
shown in blue (dashed). The circles correspond to proton induced showers and the squares have as primary particle 
iron. The error bars represent the RMS of the corresponding distribution. The points were artificially displaced for 
better visualization (proton log(£'/ey) = -0.05 and iron log(£'/ey) = +0.05. The showers were generated using 
QGSJet-II as high energy hadronic interaction model and with 6 - 40°. 



Clearly is giving information about AX'^, as it has been obtained starting from X - Xl^ax, 
with no memory of Xi . It can be measured with ground detectors only. can be combined with 
X'^-i^ax to obtain Xi similarly to what is made for the electromagnetic case fv^. 

The relation between the shapes of the electromagnetic and the muon production profile of 
each individual event can be seen in fig. |9l Here, (U-'"-,R^) are fixed to its corresponding average 
values, and most of the information is kept in a single, most sensitive variable. The correlation is 
rather strong in the most populated region of (R'^-'"-,L^). Moreover, it is almost independent of the 
primary mass composition, making it very useful for hybrid analysis. Indeed, whenever one of the 
profiles is measured accurately a prediction of the other profile shape can be established. 
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Figure 8: Shape parameters of the muon production profile for different primaries: proton (red/full) and iron 
(blue/dashed), at E = lO'^ eV, and with QGSJet-II. In (a) is shown W parameter distribution while is shown 
in (b). 



4. Dependence on the energy cutoff 

The previous studies were done with CONEX generated profiles, including a muon energy 
thresholqjof 1 GeV. The final threshold for observation depends not only on the detector used, but 
also on the zenith angle of the shower, and within the same shower will be different for different 
depths. It is worth verifying how do the previous conclusions depend on the threshold used. 

Fig. [TO]is done for proton showers generated withCORSIKA [1] (version 6.980), which allows 
the setting of lower energy thresholds (but is much slower, so there is less statistics). This version 
was modified in order to obtain the muons at production. We first checked that the results obtained 
with CORSIKA and CONEX are compatible at 1 GeV. From Fig. |10(a)|is clear that N'^^^x changes 



drastically with the cut considered, and also X'l^^ax has some variation, since low energy muons can 
be produced until much later. However all the profiles are still described by a Universal Shower 
Profile, i.e. using eq. [B which can still be described by and R^'. From Figure p^(b)| we can infer 
that gets smaller as the energy increases, and the gets slightly larger. 

To extract the full information about the muon production profile, there will be need to invert 
the propagation effects, taking into account the muon energy spectra and transverse momentum d^. 
It will also be useful to have a detector which can select muons with a minimum energy threshold. 
The systematic uncertainties coming from the translation between the modeled distribution and the 
detected one will thus be minimized. These systematic uncertainties can be estimated directly with 
the data by comparing different detection conditions, selected by zenith angle bins or by individual 



^minimum value permitted in CONEX simulations. 
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Figure 9: Correlation between the electromagnetic profile shape (characterized by R'^'" ) and the muon production 
profile shape (represented by V). The results are shown for proton (red/full) and iron (blue/dashed) induced showers, 
at E = 10'^ eV, and with QGSJet-II. The shaded area shows the one region with one sigma. 
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detectors seeing the same event. In any case, the ideas of universality and the new composition 
variables are still valid, even if the method has to be calibrated for different detection conditions 
using data. 




Figure 10: Average muon production profile as a function of the muon energy threshold, for proton induced showers at 
E = 10''* eV, with QGSJet-Il. The CORSIKA showers ai-e the average of 100 events while CONEX has 50000 events. 
In (a) the profiles are displaced to the maximum depth of the profile with E^^^ > 0. 1 GeV. In (b) average profiles are 
shown in (X',N') coordinates. 



5. Total Number of Muons and Muons from the First Interaction Point 

The integral of a Gaisser-Hillas profile is readily obtained from J N'^dX = VStt • N,„ax • L'' ■ 
f{R^), in which f{R^) is a small correction of around 1.02 for ~ 0.5. This gives the energy in 
electromagnetic profiles, and the total number of produced muons (above a given energy threshold) 
for the muon production profiles. So, by using the above formalism, the total number of muons 
can be obtained with a small uncertainty, just by detecting the region of shower maximum, in an 
event-by-event way. Moreover, as seen before, at first order, W is constant, and consequently N'^^ax 
can be used itself as a measurement of the total number of muons produced during the shower 
development. 

This is reasonably diff"erent from the usual muon counting methods, which deal with the muons 
observed at ground, more dependent on detection conditions and harder to relate to the underlying 
physics. Of course, for this the muon transport model has to take into account the produced muon 
spectrum, keep track of the muon decays through the atmosphere and detection efficiencies have 
to be considered. 

On the other hand, when working in X' , the shape variables L and R become a measurement of 
AX; together with X^ax they can be used to calculate the point of first interaction, X^ = X^ax - AX. 
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For the muon production profile, R'^ can be fixed, and W can be measured close to the shower 
maximum, so that at least a part of the variation on AX can be accounted for most of the measured 
events. In addition, the MPD distribution starts much earlier than the energy deposit profile with 
muons being produced directly from the first few interactions and with a much steeper rise. 

The fraction of produced muons that will decay before reaching the ground level increases 
with the zenith angle 6. Hence, the energy threshold for the muon to reach the ground depends 
also on this quantity, 6. Studies of the distribution of the first muons as a function of 6 can thus 
give some insight about the parent pion energy distributions in the high energy interactions. 




Figure 1 1 : Electromagnetic (red/full) and muon production (blue/dashed) longitudinal profiles of a single typical 
proton induced EAS at E = 10'"^ eV. The shower was produced using QGSJet-II as high energy hadronic interaction 
model. Both profiles were normalized to the area. The first interaction point is shown as a red star. 



6. Conclusions and Prospects 

The muon production longitudinal profile of air showers is a true characteristic of the shower 
(independent of detection conditions), described by the same parametrization as the electromag- 
netic profile. It gives new primary mass composition variables which are fairly independent of the 
high energy hadronic interaction model: the X'^l^ax and the shape variable L^. These variables can 
be combined to obtain the point of first interaction. The normalization of the profile gives also 
access to the total number of produced muons, which is known to be an important variable, both 
for primary composition and high energy hadronic interaction model studies. 
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Thus the muon production profile enclosures important information about the primary com- 
position provided that the experimental systematics are under control, as well as muon energy 
spectrum, transverse momentum and other elements that play a role on the muon transport. 

Joining all the information with the electromagnetic profile will give rise to extra variables, 
mostly sensitive to the shower development characteristics, and to a more precise understanding 
in terms of energy distribution along the shower development. With the higher number of available 
observables, using profiles which are independent of the detection conditions and more directly 
related to the hadronic cascade, cosmic rays become increasingly useful for the study of particle 
physics at the highest energies. 
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